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ABSTRACT 

We report the detection of diffuse O VI AA1032, 1038 emission in a 29-ksec observation centered on 
the Coma Cluster (I = 57.6, 6 = +88.0) and an 11-ksec observation toward Virgo (/ = 284.2, 6 = +74.5) 
through the 30" x 30" aperture of the Far Ultraviolet Spectroscopic Explorer (FUSE). The emission lines 
have a redshift near zero and are thus produced by gas in our own Galaxy. Observed surface brightnesses 
are 2000 ± 600 photons cm~ 2 s _1 sr _1 for each of the O VI components in the Coma spectrum, and 
2900+700 and 1700+700 photons cm" 2 s' 1 sr" 1 for the 1032 and 1038 A lines, respectively, toward Virgo. 
These features are similar in strength to those recently observed in the southern Galactic hemisphere 
(I = 315.0, b = -41.3) in an ~ 200 ksec FUSE observation. From a FUSE spectrum of M87, we find that 
N(0 VI) toward Virgo is (1.4 ± 0.8) x 10 14 cm' 2 . By combining emission- and absorption-line data for 
this sight line, we estimate the physical parameters of the emitting gas. 

Subject headings: Galaxy: general — ISM: general — ultraviolet: ISM 



1. INTRODUCTION 

Emission spectroscopy at extreme- and far-ultraviolet 
wavelengths will provide the key to our eventual under- 
standing of the origin, distribution, and physical processes 
of the hot interstellar medium (ISM). Observing diffuse 
interstellar emission at these wavelengths has, however, 
proven technically challenging, and only a handful of de- 
tections have been obtained to date (Martin & Bowyer 
1990; Dixon, Davidson, & Ferguson 1996). Thus, the re- 
cent observation of diffuse Galactic O VI emission with the 
Far Ultraviolet Spectroscopic Explorer (FUSE) by Shcl- 
ton et al. (2001) is both exciting, because of what it can 
tell us about the state and distribution of hot gas in the 
Galaxy, and frustrating, because of the long integration 
times (~ 200 ksec) required. We are therefore pleased to 
report the serendipitous detection of diffuse Galactic O VI 
emission in two (relatively) short FUSE observations of 
the Coma and Virgo Clusters. Combining our emission- 
line intensities toward Virgo with the O VI column density 
measured toward M87, we constrain the physical proper- 
ties of the emitting gas. These results bode well for future 
efforts to map the sky at far-UV (FUV) wavelengths. 

2. OBSERVATIONS AND DATA REDUCTION 

FUSE comprises four separate optical systems. Two em- 
ploy LiF optical coatings and are sensitive to wavelengths 
from 990 to 1187 A, while the other two use SiC coatings, 
which provide reflectivity to wavelengths as short as 905 A. 
The four channels overlap between 990 and 1070 A. For a 
complete description of FUSE, see Moos et al. (2000) and 
Sahnow et al. (2000). 

1 Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE is operated for NASA by the Johns 
Hopkins University under NASA contract NAS5-32985. 

2 Current address: Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218 
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The FUSE spectrum of the Coma Cluster was obtained 
in 17 separate exposures on 2000 June 18 and 19. Each 
exposure was centered on 12 h 59 m 49!0, +27°57'46" (J2000, 
or I — 57.61,6 = +87.96 in Galactic coordinates), near 
the center of the Coma Cluster. The total exposure 
time was 28608 s, with 23553 s obtained during orbital 
night. We use the entire 29-ksec data set in our anal- 
ysis. Our spectrum of the Virgo Cluster is a combi- 
nation of data from two locations near the cluster cen- 
ter. Two exposures, centered on 12 h 31 m 07!3, +12°23'46" 
(I = 284.03,6 = +74.52) and totaling 2242 s, were ob- 
tained on 2000 June 13. Seven exposures, centered on 
12 h 31 ro 13!4, +12°22'10" (I = 284.15,6 = +74.50) and to- 
taling 8688 s, were obtained on 2000 June 17. The total 
integration time is 10,930 s, all of it during orbital night. 
All observations were made through the 30" x 30" low- 
resolution (LWRS) aperture. The data-reduction proce- 
dures applied to these data are discussed in detail in Dixon 
et al. (2001) and are not repeated here. 

While the resolution of the FUSE spectrograph is ap- 
proximately 15 km s _1 for a point source, diffuse emission 
filling the LWRS aperture yields a line profile that is well 
approximated by a top-hat function with a width of ~ 106 
km s . Our spectra are quite faint, so we bin the Coma 
data by 8 detector pixels, or about 15.5 km s . Approxi- 
mately 7 of these bins are required to span a diffuse emis- 
sion feature, providing ample sensitivity to its line profile. 
We bin the data from the shorter Virgo observation by 
16 pixels. We distinguish between detector pixels and 8- 
or 16- pixel bins throughout this Letter. Because of the 
wide disparity in the effective areas of the various detector 
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Table 1 
O VI Emission Features 





Surface 


Intrinsic 


LSR 


Wavelength 


Brightness 


FWHM 


Velocity 


(A) 


(ph cm~ 2 s~ 1 sr _1 ) 


(km s _1 ) 


(km s _1 ) 


Coma 


1031.83 ±0.03 


2000 ± 600 


23 ± 55 


-17 ± 9 


1037.70 ± 0.06 


2000 ± 600 


75 ±60 


±35 ± 17 


Virgo 


1032.11 ±0.05 


2900 ± 700 


< 80 


±84 ± 15 


1037.55 ± 0.06 


1700 ± 700 


< 110 


±12 ± 17 



1030 



1035 
Wavelength (A) 



1040 



Fig. 1. — FUSE spectra of the (top) Coma and (middle) Virgo 
Clusters showing diffuse Galactic O VI AA1032, 1038 emission. The 
Virgo spectrum also shows Galactic C II* A1037.02 emission. Ter- 
restrial airglow features are marked with ©. The Coma spectrum is 
binned by 8 detector pixels, the Virgo spectrum by 16. The data, 
from detector LiF 1A, arc presented as histograms and arc over- 
plotted by our best-fit models. The observed continuum level is 
consistent with the dark-count rate determined from unilluminated 
regions of the detector. The round shape of the model emission fea- 
tures in the Virgo spectrum (which are in fact top-hat functions) is 
an artifact of our sparse sampling of the model lines. Bottom: FUSE 
spectrum of M87, a weighted mean of the LiF 1A and LiF 2B spec- 
tra, binned by 16 detector pixels. The error spectrum is overplotted. 
Heavy lines are the best-fit O VI absorption model, with N(0 VI) 



1.4 X 10 14 cm" 



and a Gaussian fit to the O I A1039.23 line, 



which wc adopt as the absorption profile for O VI. 

segments (Sahnow et al. 2000), combining faint spectra 
from different segments docs not significantly improve 
their signal-to-noise ratio. We thus use data only from the 
segment with the highest effective area at the wavelength 
of interest. 

The FUSE flux calibration, based on theoretical models 
of white-dwarf stellar atmospheres, is believed accurate to 
about 10% (Sahnow et al. 2000). Corrections to the nomi- 
nal FUSE wavelength scale are derived from the measured 
positions of airglow features in the day-time Coma spec- 
trum and are good to about 0.01 A. Error bars are assigned 
to the data assuming Gaussian statistics, then smoothed 
by 9 bins to remove small-scale features in the error spec- 
trum without significantly changing its shape. Segments 
of the flux- and wavelength-calibrated LiF 1A spectra of 
each cluster, showing the region about O VI AA1032, 1038, 
are presented in Fig. 1. Upper limits on the rcdshiftcd 
emission from both clusters are presented in Dixon et al. 
(2001). 

3. SPECTRAL ANALYSIS 

In our analysis, model spectra are fit to the flux- 
calibrated data using the nonlinear curve-fitting program 

3 The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy Observatories, which is operated by 
the Association of Universities for Research in Astronomy, Inc., (AURA) under cooperative agreement with the National Science Foundation. 



Note. — Error bars are statistical and represent a change in \ 2 of 1.0. 
Uncertainties in the area of the LWRS aperture and the FUSE flux 
calibration contribute an additional 14% systematic uncertainty to 
the quoted line fluxes. For lines narrower than our top-hat function, 
1-cr upper limits to the intrinsic FWHM are quoted. 



SPECFIT (Kriss 1994), which runs in the IRAF 3 environ- 
ment, to perform a x 2 minimization of the model param- 
eters. The observed profile of a diffuse emission feature 
represents a convolution of its intrinsic profile, assumed 
to be a Gaussian, with the 106 km s _1 top- hat function 
discussed above. Free parameters in the fit are the level 
and slope of the continuum (assumed linear) and the in- 
tensity, wavelength, and intrinsic FWHM of each emission 
line. The observed continuum level is consistent with the 
dark-count rate determined from unilluminated regions of 
the detector (Dixon ct al. 2001). 

The model that best fits our Coma spectrum between 
1030 and 1040 A yields a x 2 of 166.7 for 173 degrees of free- 
dom, including the O I A1039.23 airglow line. The Virgo 
spectrum, taken at night, does not show the O I feature, 
but does show emission from C II* A1037.02. (Signifi- 
cant C II A1036.34 emission is not expected; Shelton et al. 
2001). The best-fit model returns a x 2 of 80.4 for 85 de- 
grees of freedom. These models are overplotted on the 
data in Fig. 1. Parameters of the best-fit O VI lines are 
presented in Table 1, those of the C II* line in Table 2. 
In both tables, error bars for each model parameter are 
determined by increasing the best-fit value of that param- 
eter, while re-optimizing the other model parameters, until 
X 2 increases by 1.0 (corresponding to a l-er deviation for 
one interesting parameter; Avni 1976). Because our wave- 
length scale is based on the measured positions of airglow 
lines, derived velocities are geocentric, and conversion to 
the Local Standard of Rest (LSR; Mihalas & Binney 1981) 
is straightforward. 

We use the arguments of Shelton et al. (2001) to reject 
four alternative explanations for the emission observed in 
the 1030-1040 A region of the LiF 1A spectrum: scattered 
light in the cross-dispersion direction on the detector, at- 
mospheric airglow emission, solar O VI contamination, and 
H2 fluorescence. Furthermore, none of the observed fea- 
tures can be identified as rcdshiftcd emission from either 
background cluster. We conclude that the observed fea- 
tures represent diffuse Galactic emission. 

To search for statistically-significant emission at other 
wavelengths, we bin each spectrum to the instrument res- 
olution and identify individual resolution elements with 
a signal-to-noise ratio greater than 3. Each such peak 
corresponds to a single emission feature. Returning to 
our nominally-binned spectrum, we fit each feature with 
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Table 2 

Additional Spectral Features 



Wavelength 


Surface Brightness 


Intrinsic FWHM 


Detector 




(A) 


(ph cm~ 2 s _1 sr _1 ) 


(km s- 1 ) 


Segment 


Identification 


Coma 


977.04 ± 0.08 


3500 ± 2400 


43 ± 140 


SiC 2A 


C III A977.020 


1014.92 ± 0.09 


2600 ± 700 


120 ± 50 


LiF 1A 


S III A1015.505 


1122.17 ±0.04 


1600 ± 500 


< 40 


LiF 2A 


C I A1122.260 


1138.96 ±0.02 


1800 ± 500 


< 15 


LiF 2A 


Artifact? 


1144.94 ± 0.06 


2000 ± 600 


< 50 


LiF 2A 


Fe II A 1144. 938 


Virgo 


975.17 ±0.04 


6400 ± 1500 


< 40 


SiC 2A 




977.15 ±0.12 


7700 ± 3300 


117 ± 130 


SiC 2A 


C III A977.020 


1037.04 ± 0.07 


1700 ± 700 


< 110 


LiF 1A 


C II* A1037.018 



Note. — See note to Table 1. To minimize contamination by solar emission, upper limits for 
the SiC 2A channel are derived only from data obtained during orbital night. 



a linear continuum and model emission line. Because our 
data are rather noisy, we employ two measures of a fea- 
ture's statistical significance. First, we determine the 1-er 
uncertainty in its measured flux using the A\ 2 = I test 
described above; second, we test whether deleting the fea- 
ture from our best-fit model raises x 2 by at least 9 (corre- 
sponding to a 3-cr deviation for one interesting parameter, 
in this case the flux in the line; Avni 1976). We list in 
Table 2 those features that are significant at the 3-ct level 
by both tests, with three exceptions: First, we discard any 
features detected in known stray-light stripes (see Dixon 
et al. 2001). Second, we include the C II* A1037 line in 
the Virgo spectrum to facilitate comparison with a similar 
feature seen by Shelton et al. (2001). Third, we include 
the C III A977 feature, not seen by Shelton et al., because 
it is present in both of our spectra. 

Two lines presented in Table 2 remain unidentified. Be- 
cause it is so narrow, we suggest that the Al 138.96 feature 
in the Coma spectrum may represent a statistical fluctua- 
tion or detector artifact. The A975.17 feature in the Virgo 
spectrum might be redshifted Ly 7, but its wavelength 
is too short: the line would appear at 976.04 A at the 
redshift of Virgo (z = 0.0036; Ebeling et al. 1998) and 
976.78 A at that of M87 (z = 0.00436; Smith et al. 2000). 
Furthermore, most excitation mechanisms would produce 
more Ly (3 than Ly 7 emission, and redshifted Ly f3 is not 
seen in this spectrum (Dixon et al. 2001). 

4. n(o vi) towards m87 

Given the emission intensity and column density of an 
emitting species along a given line of sight, one can derive 
the electron density and thermal pressure of the emitting 
plasma (Shull & Slavin 1994). We have the O VI intensity 
toward Virgo; to constrain its column density, we consider 
a FUSE spectrum of M87, the central galaxy of the Virgo 
Cluster. 

Two night-time exposures, totaling 3611 s, were ob- 
tained on 2000 June 13 with the LWRS aperture centered 
on 12 h 30 m 49 i !l, +12°23'31" (I = 283.77,6 = +74.49), 
about 7 arcsec from the galaxy's core. To reduce the data, 
we employ the CALFUSE pipeline (Sahnow et al. 2000; 
Oegerle, Murphy, & Kriss 2000), version 1.8.7, making 
two changes to the default parameters: first, we apply the 
pulse-height cut used for our cluster data, excluding pho- 
ton events with pulse heights less than 4 or greater than 15 
in standard arbitrary units; second, we rescale the assumed 
background count rate (which depends on the pulse-height 



cut) until the residual continuum in the high- and medium- 
resolution apertures (which are essentially blank fields) is 
no more than 1-2% of the continuum in the LWRS aper- 
ture. To maximize the signal-to-noise ratio, the LiF 1A 
and LiF 2B spectra are shifted to a common wavelength 
scale (using the correction derived for Coma; see Dixon 
et al. 2001), weighted by their respective effective areas, 
averaged, and binned by 16 detector pixels. The resulting 
spectrum is presented in Fig. 1. The strongest absorp- 
tion features are due to Galactic C II A1036.34 and O I 
A1039.23. 

We generate a set of synthetic O VI absorption-line mod- 
els using an ISM line- fitting package written at U.C. Berke- 
ley by M. Hurwitz and V. Saba. Wavelengths, oscillator 
strengths, and other atomic data are taken from Morton 
(1991). Given the column density, Dopplcr broadening 
parameter, and a line list for each component, the pro- 
gram computes a Voigt profile for each absorption feature 
and outputs a high-resolution spectrum of r versus wave- 
length. We adopt the Dopplcr parameter b = 23 km s _1 
appropriate for a 5 x 10 5 K gas. The model spectra are 
convolved with a Gaussian of FWHM = 80 km s _1 (the 
observed width of the O I A1039.23 feature) and rebinncd 
to the resolution of our M87 spectrum. 

Using SPECFIT, we fit our synthetic absorption profiles 
and a linear continuum to the 1030-1034 A region of our 
M87 spectrum. The wavelength of the absorption fea- 
ture is fixed at 1032.11 A, the observed wavelength of the 
O VI A1032 emission in the Virgo spectrum. The best-fit 
model (plotted in Fig. 1) corresponds to a column density 
N(0 VI) = (1.4 ± 0.8) x 10 14 cur 2 , where the error bars 
represent A% 2 = 1. The feature is statistically significant 
at the 2-cr level, but appears both weaker and narrower 
than the C II and O I lines. Our quoted column density 
may best be considered an upper limit. 

Correcting the observed O VI doublet intensity I a ~ 
5000 photons cm~ 2 s _1 sr _1 for scattering due to dust (a 
factor of about 1.5; see below) and self-absorption within 
the emitting cloud (a factor of about 2, though r differs 
for the two components; see Shelton et al. 2001), yields an 
intrinsic O VI intensity toward Virgo of U ~ 15, 000 pho- 
tons cm~ 2 s _1 sr _1 . Using equation (5) of Shull & Slavin 
(1994), we combine h with N(0 VI) towards M87 to derive 
the electron density n e of the emitting plasma. Over the 
temperature range for which O VI is an important plasma 
constituent, 5.3 < logT < 5.8, n e is nearly constant, vary- 
ing between 0.046 and 0.048 cm -3 (though not monotoni- 
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cally), while the thermal pressure (Pth/k — l.92n e T) rises 
from 20,000 to 59,000 K cm~ 3 . This calculation assumes 
that the O Vl-emitting plasma has a constant emissivity 
per ion and that scattering by non-emitting species is neg- 
ligible. If N(0 VI) is considered an upper limit, then n e 
and Pth/k are higher than the derived values. 

5. DISCUSSION 

Shelton et al. (2001) detect O VI emission at 1032 and 
1038 A with intensities of 2930 ± 290 and 1790 ± 260 pho- 
tons cm~ 2 s _1 sr _1 , respectively, in an ~ 200 kscc FUSE 
observation centered on I = 315.0, & = —41.3. Though 
nearly antipodal, our sight lines show similar O VI inten- 
sities, suggesting that diffuse O VI emission with an ob- 
served doublet intensity of ~ 5000 photons cm~ 2 s _1 sr _1 
may be a general feature of high-Galactic-latitude sight 
lines. 

Both the Coma and Virgo Clusters show up clearly 
in the ROSAT |-kcV surface-brightness map of Snowden 
et al. (1997). The plasma responsible for the observed 
soft -X-ray (SXR) flux is likely to be somewhat hotter than 
that producing our O VI emission. Also prominent in this 
map is the North Polar Spur (NPS), thought to be gas 
shock-heated to 10 6 K by stellar winds and supernovae 
from hot stars in the Sco-Cen association (Weaver 1979; 
Heiles 1984). The Virgo Cluster lies near the edge of the 
SXR emission associated with the NPS. 

The infrared emission maps of Schlegcl, Finkbeiner, & 
Davis (1998) show that the Coma Cluster lies in a region 
of the sky nearly free of dust, with a reddening of only 
E(B — V) = 0.008. Assuming the extinction parameteriza- 
tion of Cardclli, Clayton, & Mathis (1989) and R v = 3.1, 
this reddening corresponds to an attenuation of ~ 10% at 
1035 A. The infrared sky is more complex in the direction 
of the Virgo Cluster, as the cluster lies on the edge of a 
broad dust lane that may be associated with the NPS. The 
extinction toward Virgo is E(B-V) = 0.030. Any O VI 
AA1032, 1038 emission originating beyond this obscuring 
material will be reduced in intensity by ~ 34% (hence our 
reddening correction of 1.5). 

The column density N(0 VI) derived from our M87 
spectrum falls at the low end of the range N(0 VI) 



- (1.4-5.0) x 10 14 cm" 2 obtained by Savage et al. (2000) 
for sight lines through the Galactic halo. Such a low value 
of N(0 VI) suggests that our Virgo sight line does not in- 
tercept significant O Vl-emitting plasma associated with 
the North Polar Spur. If so, then parameters derived 
for this sight line may be typical of the Galactic halo. 
Given that the total pressure (thermal + nonthermal) in 
the Galactic plane may be as low as ~ 25,000 K cm -3 
(McKee 1993), the thermal pressure we derive for the halo 
seems quite high. The derived thermal pressure would 
rise if N(0 VI) is lower than our best-fit value, but would 
fall if the cirrus clouds responsible for dust extinction lie 
beyond the O Vl-emitting plasma. Additional FUSE ob- 
servations of M87, already scheduled for Cycle 2, should 
yield a greatly-improved determination of N(0 VI). 

As seen in Table 1, the LSR velocities of the O VI 1032 
and 1038 A emission features differ significantly from one 
another in both the Coma and Virgo spectra. Neither in- 
strumental calibration uncertainties nor contamination by 
emission from other species can explain this effect. More 
likely, the apparent velocity differences are optical-depth 
effects. Because the optical depth of the 1038 A line is 
half that of the 1032 A line, the 1038 A flux that we ob- 
serve may have been emitted in regions far beyond those 
responsible for the 1032 A flux. O VI component velocities 
from multiple sight lines may thus be used to map velocity 
fields in the Galactic halo, even when individual velocity 
components are not resolved. 

Our serendipitous detection of diffuse O VI emission 
bodes well both for continued FUSE observations of se- 
lected sight lines and for proposed missions, such as 
SPEAR (Spectroscopy of Plasma Evolution from Astro- 
physical Radiation; Edelstein, Korpela, & Dixon 2000), 
which seek to map the entire sky at FUV wavelengths. 
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